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While the modification of jets by a medium as created in ultrarelativistic heavy-ion collisions 
is often thought of as a phenomenon affecting a hard near on-shell parton and hence treated as 
partonic energy loss, a more realistic view is to think of the medium effect as a relatively small 
correction to the Quantum Chromodynamics (QCD) vacuum radiation pattern which reduces the 
initial high virtuality of hard partons and leads to the development of a shower. The uncertainty 
relation can be used to argue that a significant part of this shower is developed even before a medium 
can be formed. The initial virtuality of a shower-initiating parton is reflected in the final state in 
the measured invariant mass of a jet — high mass jets are characterized by more branchings, higher 
multiplicity and a wider angular structure. Since most of this structure is determined before the 
medium affects the jet, selecting a particular jet mass range strongly biases the partonic configuration 
that enters the medium, and thus the medium is expected to modify high mass jets more strongly 
than low mass jets. In this work, this scenario is explored using the in-medium shower evolution 
code YaJEM and a possible strategy for a measurement is suggested. 
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I. INTRODUCTION 



A hard QCD process typically results in highly virtual 
final state partons which subsequently shed virtuality in 
a characteristic radiation pattern accessible in perturba- 
tive QCD (pQCD) and evolve into parton showers which 
after hadronization become experimentally accessible as 
a jets, i.e. collimated sprays of hadrons. In a virtuality- 
ordered shower description such as e.g. the PYSHOW 
algorithm [l| , the radiation pattern is generated as a se- 
ries of 1 — > 2 splittings a — » b, c where the parton virtual- 
ities decrease Qb,Q c *C Q a - Energy-momentum balance 
then dictates that the lost parton virtuality leads to a 
transverse momentum difference between the daughter 
partons b and c, however the invariant mass of the whole 
system is conserved. The branchings terminate as soon 
as parton virtuality reach a non-perturbative scale O(l) 
GeV, at which point pQCD is no longer applicable. The 
virtuality Qi of the shower initiating parton thus gov- 
erns both the average number of steps in the branching 
chain and the transverse momentum distribution of the 
jet fragments. 

Experimentally jets are defined by applying a cluster- 
ing algorithm to the final state hadrons, which typically 
involves (explicitly or implicitly) an angular distance R 
of hadrons with respect to the jet axis and a minimum 
constituent transverse momentum Pt cut as well as pos- 
sible PID effects (the latter two reflecting the fact that 
not all hadron species reach the calorimeter or leave a 
signal). This means that the strict relation between jet 
invariant mass Mj et and initial virtuality QiS broken — 
hadrons at large angles <fi > R with the jet axis carry 
the imprint of the original parton virtuality but are not 
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counted as part of the jet and hence do not contribute to 
jet mass. This implies that any realistic calculation can 
not use Qi as a theoretical proxy for Mj e t but must take 
into account the clustering procedure. 

The uncertainty relation allows to estimate the 
timescale at which a splitting takes place as E/Q 2 with 
E the intermediate parton energy scale. Since the ini- 
tial virtuality is of order of the hard momentum scale 
Qi ~ 0(p T ) ~ 0(100 GeV) for LHC kinematics, the 
timescale associated with the first splittings is paramet- 
rically ~ 0(0.01 fm)), i.e. far smaller than the medium 
formation time tq ~ 0(0.2 — 0.5 fm). What is potentially 
modified by the medium when the hard process takes 
place in a heavy-ion collision is thus the parton configu- 
ration as it has been created by the initial branchings by 
the time the medium is formed, and as discussed above, 
the nature of this configuration is correlated with Qi or 
M jet . 

The importance of going beyond an energy loss ap- 
proximation for a leading on-shell parton and comput- 
ing this in-medium shower evolution has been pointed 
out already in Intuitively one might expect that the 
medium modification of a shower depends on the mul- 
tiplicity of showering partons by the time the medium 
is forming as each of the shower partons may interact 
with the medium, and thus the total momentum trans- 
fer between medium and shower partons scales with the 
number of available partons. As a result, one may ex- 
pect that jets with large Mj et show on average more sub- 
stantial medium modifications than jets with small Mj et . 
In this work, we explore this idea using the in-medium 
shower evolution Monte Carlo (MC) code YaJEM 0-H] 
in its latest version YaJEM-DE [6J, test if the effect is 
robust against the bias introduced by clustering and sug- 
gest possible experimental strategies to measure it. 
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II. THE MODEL 



In order to illustrate the consequences of the jet mass 
effect most directly, let us consider the simple case of 
a high pt quark with fixed energy undergoing shower 
evolution and fragmentation in the vacuum and in the 
presence of a medium. Experimentally, such a situation 
can be realized approximately in 7-I1 correlations where 
a trigger condition on the 7 restricts the momentum dis- 
tribution of the recoiling parton (predominantly a quark) 
to a comparatively narrow region. We are interested in 
studying properties of the jets resulting from the QCD 
evolution of these quarks in vacuum and in medium. 

In order to simulate the observable effects of Qi , back- 
to-back quark pairs at given energy E = 200 GeV are 
generated as the final state of a hard process with the 
virtuality distribution of the quarks given by the initial 
splitting of the hard transient state into the final state 
of a back-to-back quark pair. One of the quarks is then 
selected and evolved with YaJEM down in virtuality to 
the hadron level. The resulting events are clustered using 
the anti-fey algorithm from the Fast Jet package Q and 
binned according to Mj. We then study the longitudinal 
momentum distributions of hadrons inside the clustered 
jets in small and large Mj bins and compute the ra- 
tio of the distributions which corresponds to Iaa(Pt) in 
the experimental realization as 7-h correlation. In addi- 
tion, we also study the transverse distribution of hadrons 
dN/dcf>(PT) as a function of the angle with the jet axis cf> 
in two different Pt bins. The experimental observation 
would involve a 7-h coincidence using the jet axis that 
has been determined in the previous step (the need to 
do a correlation arises since the jet cone size R otherwise 
restricts the angular region which can be probed). 

The core of the calculation, i.e. the evolution of par- 
ton showers in the medium is computed with the MC 
code YaJEM, which is based on the PYSHOW code [J 
which in turn is part of PYTHIA [8j. It simulates the 
evolution from an initial parton with virtuality Qi to a 
shower of partons at lower virtuality in the presence of 
a medium. In the absence of a medium, YaJEM by con- 
struction reproduces the results of PYSHOW. A detailed 
description of the model can be found in ■ Here we 
use the version YaJEM-DE Q which is one of the best- 
tested theoretical models available for in-medium shower 
evolution and gives a fair account of a large number of 
high P T observables both at RHIC and LHC [Ml. 

In YaJEM-DE, the medium is characterized by two 
transport coefficients, q and e. Here, q parametrizes 
the virtuality growth of a parton per unit pathlcngth 
and leads to medium-induced radiation whereas e de- 
scribes the energy loss of propagating partons into non- 
perturbative medium modes. These transport coeffi- 
cients are assumed to be related to the energy density 
e of the hydrodynamical medium as 



q[e}(0 = K[K D ] ■ 2 ■ [e(C)] 3/4 (coshp(C) - sinhp(C) cost/;) 

with Q the position along the path of the propagat- 
ing parton, p the medium transverse flow rapidity, ip the 
angle between parton direction and medium flow vec- 
tor and K, Kjj two parameters regulating the strength 
of medium-induced radiation vs. direct energy loss into 
the medium. As in Q, the free parameters K,K]j are 
adjusted such that the energy loss into non-perturbative 
modes is a 10% contribution to the total as constrained 
by a number of other observables 0] . 

Given (q((), e(C)) along the parton path, YaJEM com- 
putes the in-medium partonic shower evolution and 
hadronizes the result using the Lund string model (l3j 
so that the fragmentation remnants of the initial hard 
parton can be analyzed on the hadron level. However, a 
scaling law identified in [4j allows to simplify the compu- 
tation by ordering shower according to their total line- 
integrated virtuality gain AQ 2 ot — J d(q((), i.e. the pre- 
cise functional form of q(() is significantly less important 
that the line-integrated value. 



III. RESULTS 

In Fig. [TJ the ratio of medium-modified over vacuum 
fragmentation functions (equivalently away side Iaa in a 
7-h correlation) is shown for a fixed AQ 2 ot = 10 GeV 2 
which represents the medium modification for an average 
path in the medium as created at the LHC (see [l4| for 
details of the medium evolution model in which this is 
derived) in a low jet mass range Mj et = 2 — 8 GeV and 
a high jet mass range Mj et = 14 — 20 GeV. 
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FIG. 1: Ratio of medium-modified over vacuum fragmenta- 
tion function for a 100 GeV quark for two different jet mass 
ranges (lines to guide the eye). Jets have been clustered with 
the anti-fcr algorithm with R — 0.4. 

(As a side remark, note that the use of the average 
AQ 2 ot is not justified in a detailed comparison with mea- 
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sured data which reflect the fluctuations around the av- 
erage in a non-trivial way, but here it is a meaningful 
procedure to establish qualitative features and approx- 
imate magnitude of the effect. This is because unlike 
other observables which involve a convolution with the 
parton spectrum, such as the nuclear modification fac- 
tor Raa, the modification of the fragmentation function 
studied here is appoximately a linear function of AQ^ ot 
and fluctuations can hence be expected to largely average 
out.) 

It is clearly visible from the figure that, as expected, 
high mass jets which map to an increased shower activity 
before the medium modification subsequently also show 
stronger medium effects whereas low mass jets remain 
fairly robust against medium modification. 
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shows the same pattern of suppression at high Pt and 
enhancement at low Pt as seen in Fig. [T] However, the 
effects are not unifrom in angle and the most dramatic 
differences between high and low jet mass and high and 
low Pt window can be observed at <j> > 0.4 (i.e. outside 
the jet finding criterion used in this study so that a dif- 
ferent experimental technique would be needed to study 
them). Again this is not unexpected as higher mass jets 
have more phase space available to populate large angles. 



IV. CONCLUSIONS 

The simple study presented here establishes that there 
is reason to expect a non-trivial interplay between the 
vacuum radiation pattern as driven by the shower- 
initiating virtuality Qi and the subsequent medium- 
induced modification of the branching pattern: in a 
virtuality-ordered shower picture, showers get modified 
proportional to the number of virtual partons present by 
the time the medium is formed. An experimental anal- 
ysis of observables differential in jet mass Mj et is thus a 
non-trivial test for models of in-medium shower evolution 
and also in principle carries the imprint of the medium 
formation time. The expected effects have been demon- 
strated to be robust against clustering (which tends to 
blur the relation between Qi and Mj et ) and they appear 
to be large enough to be observed in practice through 
e.g. 7-jet correlations. Jet mass differential observables 
offer thus an exciting future opportunity for hard probes 
studies at RHIC and LHC. 



FIG. 2: Ratio of medium-modified over vacuum angular 
hadron distribution of a 100 GeV quark for two different jet 
mass and hadron Pt ranges. Jets have been clustered with 
the anti-fcr algorithm with R = 0.4. 

A more differential look at ratios of the transverse 
medium and vacuum hadron distribution dN / d<f> with re- 
spect to the jet axis in different Pt bins as shown in Fig. [5] 
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